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Abstract 

We derive simple new expressions, in various dimensions, for the functional determinant of a 
radially separable partial differential operator, thereby generalizing the one-dimensional result of 
Gel'fand and Yaglom to higher dimensions. We use the zeta function formalism, and the results 
agree with what one would obtain using the angular momentum cutoff method based on radial 
WKB. The final expression is numerically equal to an alternative expression derived in a Feynman 
diagrammatic approach, but is considerably simpler. 
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I. INTRODUCTION AND RESULTS 



Determinants of differential operators occur naturally in many applications in mathemat- 
ical and theoretical physics, and also have inherent mathematical interest since they encode 
certain spectral properties of differential operators. Physically, such determinants arise in 
semiclassical and one-loop approximations in quantum mechanics and quantum field theory 

013 EL 

Determinants of free Laplacians and free Dirac operators have been extensively 
studied 0,0,0,0,0,0,13], but much less is known about operators involving an arbitrary 
potential function. For ordinary (i.e., one dimensional) differential operators, a general the- 
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151 ] . In this 



ory has been developed for determinants of such operators 
paper we extend these results to a broad class of separable partial differential operators. 
The result for four dimensions was found previously in jlil us ing radial WKB and an angu- 
lar momentum cut-off regularization and renormalization 171. Here we re-derive this result 
using to Zet a taction approach t0 det _ ts Qflflff and g ene r a liz e to other dnnen- 
sions. We are motivated by applications in quantum field theory, so we concentrate on the 
dimensions d = 2,3,4, but the mathematical extension to higher dimensions is immediate. 
We also compare with another expression for the four dimensional determinant , derived 
using a Feynman diagrammatic approach. 

Consider the radially separable partial differential operators 



M = -A + V{r) 



;i.2) 



where A is the Laplace operator in jR d , and V(r) is a radial potential vanishing at infinity 



as r 2 e for d = 2 and d — 3, and as r 4 e for d = 4. For d = 1, wit 



i Dirichlet boundary 



conditions on the interval [0, oo), the results of Gel'fand and Yaglom [h] lead to the following 
simple expression for the determinant ratio: 

det [M + m 2 } ip(oo) 



det [M bee + m? 



free/ 



OO 



;i.3) 



where [Ai + m 2 ] ip = 0, with initial value boundary conditions: ^>(0) = and ip'(0) = 1. 



free 



m 



The 



The function -0 is defined similarly in terms of the free operator: [M. 
squared mass, m 2 , is important for physical applications, and plays the mathematical role 
of a spectral parameter. The result (jl.3j) is geometrically interesting, in addition to being 
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computationally simple, as it means that the determinant is determined simply by the 
boundary values of the solutions of [M + m 2 ] ip = 0, and no detailed information is needed 
concerning the actual spectrum of eigenvalues. 

Now consider dimensions greater than one (as mentioned, we are most interested in 
d = 2,3,4; but the extension to higher dimensions is immediate). Since the potential is 
radial, V = V(r), we can express the eigenfunctions of Ai as linear combinations of basis 
functions of the form: 



iMO*(q(0) 



1.4) 



r (d-l)/2 

where Ym(8) is a hyperspherical harmonic Q], labeled in part by a non-negative integer I, 
and the radial function ipm (r) is an eigenfunction of the Schrodinger-like radial operator 



M (l) = - 



(* + ¥) ('+¥) 



dr 2 



+ V(r) 
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-M|^ e is defined similarly, with the potential omitted: V = 0. In dimension d > 2, the radial 
eigenfunctions have degeneracy given by 

deg(/;d) = • (1.6) 

Formally, for the separable operators in (jl.lj) — (|1.2j) . the logarithm of the determinant ratio 
can be written as a sum over I (weighted with the degeneracy factor) of the logarithm of 
one-dimensional determinant ratios, 



In 



det[.M 



m 



det [M iree + m? 



^deg(/;rf) In 



1=0 
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Each term in the sum can be computed using the Gel'fand-Yaglom result (|1.3|) . However, the 
I sum in fll.7|) is divergent, as noted by Forman Q| for the free Laplace operator in a two- 
dimensional disc. In this paper we show how to define a finite and renormalized determinant 
ratio for the radially separable partial differential operators (jl.lj) - (jl.2j) . Specifically, we 
derive the following simple expressions, which generalize (jl.3j) to higher dimensions: 



In 



det [M + m 2 
det\M~ 



free 



77T 



hi 



^(0)( 



OG 



'.free / 



d=2 




^2 m 



oo 



J°° drrV{r) 



i=i 



free / 

(0 1 



oo 
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+ J drr V [hi ^ 
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det [M + m 2 
det[.M free + m 21 

det[A4 + m 21 
det[M ire 




J °° dr rV{r) 
2(1+1) 



;i.9) 



J °° r V(r) J °° dr r 3 V^ (V + 2m 2 



2(7 + 1] 
drr 3 1/(^ + 2m 2 ^ 



In 



//■/■ 

y 



8(Z+1) 3 
+ 7 + 1 
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Here 7 is Euler's constant, and \i is a renormalization scale (defined in the next section), 
which is essential for physical applications, and which arises naturally in even dimensions. A 
conventional renormalization choice is to take \i — m in (jl.8j) - (jl.l()j) . In each of (jl.8j) - (jl.l()j) . 
the sum over I is convergent once the indicated subtractions are made. The function ipa)(r) 
is the solution to the radial equation 



[M { i) +m 2 ] ^(j)(r) 







r 



2+(d-l)/2 



as r 
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The function V'([) ee ( r ) is defined similarly, with the same behavior as r — > 0, in terms of the 
operator [-M|^ c + m 2 



Thus, in d dimensions, ip^ e (r) is expressed as a Bessel function: 



/ 9 \ J+d/2-l 



/ + - ) r 1 / 2 J, +d/2 -i(mr) 
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Notice that the results (|1.8)) - ()1.10p state once again that the determinant is determined 
by the boundary values of solutions of [M. + m 2 ] ip — 0, with the only additional information 
being a finite number of integrals involving the potential V(r). We also stress the compu- 
tational simplicity of (jl.8|) - ()1.10|) . as the initial value problem (|1.11 | ) is trivial to implement 
numerically. The d = 4 result ()1.10|) was found previously in |l(| using radial WKB and 
an angular momentum cutoff regularization and renormalization Q]. Here we present a 
different proof using the zeta function approach, and generalize to other dimensions. In 



'act , th e d = 2 and d = 3 results can also be derived using the radial WKB method of 
1 fillip. Furthermore, in Section ITTT1 we show how these results also agree with the Feynman 
diagrammatic approach, by showing that the d = 4 zeta function expression (jl.lOj) agrees 
precisely with a superficially different, and more complicated, d = 4 expression found in jl^ |. 



Finally, we note that the results in ([1.8|) - ([1.10|) are for a generic radial potential V(r). There 
are important physical applications where the potential V(r) is such that the operator Ai 
has negative and/or zero modes ^(|, in which case these expressions are modified slightly, 



as m 



1611 and as discussed below in the conclusions. 



II. ZETA FUNCTION FORMALISM 

The functional determinant can be defined in terms of a zeta function 
operator Ai. For dimensional reasons, we define 



ODD 



C[A4+m 2 ]//i 2 ( s ) — C[M+m?]{ s ) ~ ^ 



m 



2\-s 



for the 



(2-1) 



where the sum is over the spectrum of Ai, and fi is an arbitrary parameter with dimension 
of a mass. Physically, fi plays the role of a renormalization scale. Then the logarithm of the 
determinant is defined as |4|, |5, 



lndet[.M + m 2 ] = -([ M+m2]/)l2 (0) 

= - ln(/i 2 ) C[M+m2] (0) - C[ M+m ^ (°) • 

To compute the determinant ratio, we define the zeta function difference 

C( s ) = C[M+m 2 ](s) - C[A4 frcc +m2]( s ) • 



(2.2) 



(2.3) 



Thus we need to compute the zeta function and its derivative, each evaluated at s — 0. In 
general, the zeta function at s = is related to the heat kernel coefficient, 6^/2 (P), associated 
with the operator V li], |3| : 



Cp(0) = a d/2 {V) . 



(2.4) 



For the operator V 



-A — E, these heat kernel coefficients are [10] given in 1R by 
1 



ai(V) 



(47r)^/2 J M 



d d xE 



a 3/2 (V) = 
a 2 (V) 



^)d/2 J Md 2 



(2.5) 
(2.6) 
(2.7) 



5 



Thus, setting E = —V(r) — m 2 , and E iree = —m 2 , we find 



l r°° 
2 



f™drrV(r) , d = 2 



C(0) = < , d= 3 (2i 



1 J °° drr 3 V(V + 2m 2 ) , d = A 
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The derivative of the zeta function at s = 0, C'(0), can be evaluated using the relation 



to the Jost functions of scattering theory 



21 



for the application of these ideas to the 



Casimir effect see |23j. Consider the radial eigenvalue equation 

M(i)(f>(i),p = p 2 (j>(i), P , (2.9) 

where Mn\ is the Schrodinger-like radial operator defined in (jl.5j) . A distinguished role is 
played by the so-called regular solution, <f)(i) )P (r), which is defined to have the same behavior 
as r — > as the solution without potential: 

0(O, P ( r ) ~ Ji+(d-3)/2{pr) ■ (2.10) 
Here the spherical Bessel function ji+(d-3)/2 is defined as 

3l+(d-3)/2(z) = ^-y Ji+d/2-l(z) . 

The a sym ptotic behavior of the regular solution, 0«),p(r), as r — > oo defines the Jost function, 



fi(p), m 



% r a a 

0(o, P (r) ~ - [/Kp)^r+( d _3)/ 2 (p r ) -/r(p)^( d _3)/ 2 (p r )J • ( 2 - n ) 



Here h l+ , d _ z y 2 {pr) and h^ d _ 3 y 2 (pr) are the Riccati-Hankel functions 

^t+(d~3)/2( Z ) = ^^~^^l+d/2-l( Z ) ' ^«+(d~3)/2( Z ) = ~^\j~2^l+d/2-xi Z ) ■ 

As is well known from scattering theory jsiL l^ . the analytic properties of the Jost 
function fi(p) strongly depend on the properties of the potential V(r). Analyticity of the 
Jost function as a function of p for > is guaranteed, if in addition to the aforementioned 
behavior as r — > oo, we impose V(r) ~ r" 2+e for r — > 0, and continuity of V(r) in < r < oo 
(except perhaps at a finite number of finite discontinuities). For us, the analytic properties 
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of the Jost function in the upper half plane will be of particular importance because they 
are related to the shifting of contours in the complex momentum plane. 

By standard contour manipulations Q|, the zeta function can be expressed in terms of 
the Jost functions as: 

c(s) = ^^ degM) d k[k 2_ m 2 ] -s ° lnfi{ik) (212) 
1 u m 

This representation is valid for Jfe > d/2, and the technical problem is the construction of 
the analytic continuation of ()2.12j) to a neighborhood about s = 0. If expression ()2.12jl were 
analytic at s — 0, then we would deduce that 

00 

c'(o) = -E de g^ ln ^ m ) • ( 2 - 13 ) 

1=0 

From the definition (j2.11|) of the Jost function, 

where ip(i){r) is defined in (ll.llj) . Thus, the regulated expression ()2.13j) coincides with 
the formal partial wave expansion (jl.7|) . using the Gel'fand-Yaglom result (jl.3j) for each 
/. However, the expansion ()2.13|) is divergent in positive integer dimensions. In the zeta 
function approach, the divergence of the formal sum in (j2.13|) is directly related to the need 
for analytic continuation of ((s) in s to a region including s = 0. From (j2.12)) . this analytic 
continuation relies on the uniform asymptotic behavior of the Jost function fi{ik). Denoting 
this behavior by f^ sym (ik), the analytic continuation is achieved by adding and subtracting 
the leading asymptotic terms of the integrand in (J2.12j) to write 



COO = CM + Cas(s) , (2.15) 



where 



and 



= ^JdegM /*^-m^pn/i(tt)-lnH*)] > ( 2 - 16 ) 
Us) = ^M^degM) [dk[e-m 2 r s -^lnfr m W . (2.17) 

71 1=0 L 



Ultimately we are interested in the analytic continuation of £(s) to s = 0. As many asymp- 
totic terms will be included in fi Sym [ik) as are necessary to make C/( s ) as given in (j2.16j) 
analytic around s = 0. On the other hand, for ( as (s) the analytic continuation to s = can 
be constructed in closed form using an explicit representation of the asymptotic behavior of 
the Jost function, derived in the next section. 

A. Asymptotics of the Jost Function 



The asymptotics of the Jost function fi(ik) follows from standard results in scattering 



theory 



2l| . The starting point is the integral equation for the regular solution 



0(O,pO) = ji +{d - 3 y 2 (pr) + J dr' g l>p (r, r') V(r') (j>(i), p {r') 

o 

with the Green's function 



2p 



h 



l+(d-3)/2(P r )K+{d-3)/2(P r ) ^t+(d-3)/2(P r )^l+(d-3)/2(P r ) 



Asymptotically for r — > oo, 

0(O,p( r ) ~ Ji+(d-3)/2(pr) + / dr' Gi, p (r, r') V{r') (j){i), P {r') . 



(2.19) 



'2.20) 



Noting that [h„ (x)]* = h^(x), for x real, this asymptotic behavior can be written as 



(0.p 



r ~ 



^pj dr ' K+ {d - m (pr')V{r>)ct> {l) , p {r') 
o 



1 + pJ dr ' 





h 



l+(d-3)/2 



(pr) 



h 



l+{d-3)/2 



(pr) 



Comparing with the definition (J2.11J) of the Jost function, we find the following integral 
equation for the Jost function: 

oo 

/,(p) = 1 + i J dr ht +M/2 (pr) V(r) {OiP (r) . (2.21) 
o 

For the zeta function (j2.12|) we need the Jost function for imaginary argument, so we rotate 
using the Bessel function properties 



Uz) = e—^JJiz) , KJz) 



Til 



iz) 



Thus (J2.21|) becomes 



fi(ik) = l+ / drrV(r)(j)(i) iik (r)K l+d/2 -i(kr) 



(2.22) 



We define a convenient short-hand for the Bessel function index, 

, d 1 
v = I H 1 . 

2 



(2.23) 



and write the partial-wave Lippmann-Schwinger integral equation for the regular solution 

as 



= I„(Ax) + / drV [I v {kr)K v {kr') - I v {kr')K v {kr)]V{r')cj) {l)>ik {r^ . (2.24) 



This Lippmann-Schwinger equation leads to an iterative expansion for fi(ik) in powers of 
the potential V(r). For dimensions d < 4, we need at most the 0(V) and 0(V 2 ) terms of 
ln/,(iA;): 

oo oo r 

]nfi(ik) = J drrV(r)K u (kr)I u (kr) - J dr rV(r)K 2 (kr) J dr' r' V(r')I 2 (kr') + 0(V 3 ) . 



(2.25) 

This iterative scheme effectively reduces the calculation of the asymptotics of the Jost func- 
tion to the known uniform asymptotics of the modified Bessel functions K v and I v j^J]. To 
the required order, for v — > oo, k — * oo, with kjv fixed, 



I u {kr)K u {kr) ~ 1- + (1 - 6t 2 + 5t 4 ) + o(\) , 
h{kr')K v {kr) 



-v(rj(k)—r)(kr' /r)) 



l + O [ - 

v 



(2.26) 



Iv (1 + (fcr/i/) 2 ) 1/4 (l + (fcr'/z/) 2 ) 1 ^ 

where t = + (kr/v) 2 , and ^(A;) = a/1 + (kr/v) 2 + ln[(Ax/i/)/(l + + (kr/vW)}. 

The r' integration in the term quadratic in V is performed by the saddle point method |10( ■ 
We therefore define In f^ sym (ik) as the 0(V) and 0(V^ 2 ) parts of this uniform asymptotic 
expansion: 



oo 

\nfi sym {ik) = ^- J dr 



r V(r) 



kr \ ^ 



1/2 



9 



16^ 3 

1 

8^ 



dr 



r V(r) 



o 

oo 



dr 



r 3 V 2 (r) 



3/2 



+ 



kr \ ^ 



n 2 







3/2 



(2.27) 



B. Computing C/(0) 



By construction, C/( s ); defined in ()2.16|) . is now well defined at s — 0, and we find 

oo 

C/(0) = -^deg(/;rf)[ln/K^)-ln/r m (^)] • (2-28) 



z=o 



This form is suitable for straightforward numerical computation, as the Jost function fi(im) 
can be computed using (jl.llj) and (j2.14j) . while hi f^ sym (im) can be computed using ()2.27j) . 



With the subtraction of In f[ 



asym 



im, 



in ()2.28|) . the I sum is now convergent. 



However, it is possible to find an even simpler expression. It turns out that the subtraction 
in (j2.28|) is an over-subtraction. To see this, expand In f^ sym (im) into its large / behavior as 
follows: 



oo oo 

In f^ sym (im) ~ — J drrV(r) — -^—^ J dr r 3 V (V + 2m 2 



H / drrV(r) 



/mr 
V v 



i -i/2 




16u 3 



dr 



r V(r) 



o 



3/2 



OS' 



drr 3 V 2 (r) 



r 


(mr \ 2 


-3/2 ^ 




1 + (-) . 











(2.29) 



The first term is O (j), and the second is O (js), while the remaining terms are all O (js)- 
In dimensions d < 4, the degeneracy factor deg(/; d) is at most quadratic in /, and so these 
last terms are finite when summed over I in (JOED- (In fact, in d = 2 and d = 3, the O (£) 
terms are also finite when summed over I.) In the next section we show that these finite 
terms cancel exactly against corresponding terms arising in the evaluation of C as {0). Thus, 
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for C'(0) = C/(0) + Cas(O) , we only actually need to subtract the leading large I terms in 
()2.29j) . rather than the full asymptotics in ([2.27)1 . 



C. Computing C«(0) 



The explicit form of the asymptotic terms in (|2.27|) provides the analytic continuation to 
s = of C as (s), as defined in ()2.17|) . The k integrals are done using 

2" 



dk \k 2 -m 2 ]- s ^- 
ak 



( kr 

1 + 



T(s + 1 )r(i - ,) (£)" 



III 



-2s 



r(n/2) 



V mr I 



~n~ ■ (2.30) 



Therefore, we find 



Cas(s) = -^deg(/;d) 



1=0 



™ drr 1+2s V(r) \ 1 F(s + ^ 



-l-2s 







+ 



+ 



1 T(s + 



v 



-3-2s 



2r( a )r(|) + 
3r(s + |) 



s+l/2 



3-2s 



i6r( s )r(|) ^^2^+3/2 8r(s)r(|) + ^ 
5 r( a + j; 



s+5/2 



-3-2s 



i6r( s )r(|) f 1 + ^ 



s+7/2 



drr s+2s V 2 (r 



ir( a + §; 



-3-2s 



sr( s )ro 



1 + (^y 



s+3/2 



(2.31) 



We now subtract sufficiently many terms inside the I sum to ensure the analytic continuation 
°f (as{s) to s — 0. The added back terms produce Riemann zeta function terms, such as 
Cr(2s + 1), whose analytic continuation is immediate. For example, in d = 4, where v = l + l, 
and deg(Z; 4) = (/ + l) 2 = z/ 2 , the first term in ()2.3H) involves the function 

,,l-2s 



^(*) = r( 8 )r(^ £ 



.2) ,=1(1 

The analytic continuation of this function to s = is 



8+1/2 



(2.32) 



Tills) 



\)r 2s 



room; 



i/=i 



/ mr \2\ s 1 / 2 / ]_\ / mr \: 



+C*(2s - 1) - ( s + - ) (mr) 2 ( R (2s + 1) 



(2.33) 
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A straightforward computation yields the derivative at s = 0: 

-1/2 



* ( o,=|>{Hv) i ) 



— mr 
2 V ' 



In 



+ 7 + 1 



1 fmr\ 2 
1 + 2vV) 



+c«(-i: 



(2.34) 



Applying this strategy to the remaining terms in ()2.31j) leads to 

"OO 



C(o) 



d=4 







drr 3 V{V + 2m 2 ) In 



drr V(r) < - 



+ 7 + 1 



i/=i 



+ 16 ^ i/ 



/mr 
1+ — 
V v 



-3/2 



1 f°° °° -i 



1 - 


/mr 




f (v 


f 


^mr\ 2 






1 - 


/mr 




H7 



-1/2 



1 / mr \ 2 

1 + 2b 



-5/2 



5 (i+r^ 



V 1/ 



-7/2' 



-3/2 



'2.35) 



Notice that the terms involving summation over v cancel exactly against identical terms in 
C/(0) from (J2.29)) . after those terms are summed over I with the d = 4 degeneracy factor 
v 2 = (I + l) 2 . Furthermore, note that the lnr term inside the integral on the first line of 
f!2.35|) is precisely of the same form as the renormalization term in (|2.8|) . so the In \x in ()2.2}) 
combines with lnr to form the dimensionless combination ln(/xr) in (jl.lOj) . 

The analogous computation in d = 3, with degeneracy factor deg(7; 3) = {21 + 1) = 2u, 
leads to 



C(o) 



d=3 



oo 

8 to (^l) 



i + 



3/2 



1 + 



5/2 



i + 

2\ ~ 3 / 2 



7/2 



(2.36) 



In this case, all terms in (J2.36)) cancel against corresponding terms in ()2.29j) . after summing 
over I with degeneracy factor 2v in d = 3. The only remaining uncancelled term in ()2.29|) is 
the first term, which is linear in V, and is the subtraction shown in (|1.9|) . This shows that 
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in dimension d = 3 we did not actually need to expand In f^ sym (ik) to 0(V 2 ) in the first 
place; the 0(V) term would have been sufficient. 

The d = 2 case is slightly different, as we need to separate the / = term from the sum. 
Here v — I, and the degeneracy factor is 1 for I = 0, and 2 for I > 1. Thus, instead of (J2.28|) 
we have 



c;(o) 



d=2 



- In /o(zm) -^2 [In fi{im) - In ff sym {im)} 



i=i 



And in two dimensions ()2.17|) is 



(2.37) 



Cas(s) 



^^f> Jdk [k 2 -m 2 }- s -^\nfr m i 



sm(7rs 

d=2 



ik) 



Then the analogous computation in d = 2 leads to 



CL(o) 



d=2 



drr V(r) |ln ( - J +7 



f 00 1 

drr V(r) < - 



1=1 



/mr\ 2 
W 



!=1 
-3/2 



fmr 

1 + (t 



-1/2 



- 1 



/mr\- 

+ (-) 



-5/2 



:2.38) 



5 1 



-7/2' 



■if*"M£H 1+ (r)') 



-3/2 



(2.39) 



As in the d = 4 case, all terms involving / summation cancel exactly against identical terms 
in (|2.29j) . after summing those over /, with the d = 2 degeneracy factors. As in d = 3, the 
only remaining uncancelled term in (|2.29J) is the first term, which is linear in V, and is the 
subtraction shown in (jl.Sjl . This shows that also in dimension d = 2, we did not need to 
expand In f^ sym (ik) to 0(V 2 ) in the first place; the 0(V) term would have been sufficient. 



III. COMPARISON WITH FEYNMAN DIAGRAM APPROACH 

In this section we show that our zeta function computation is equivalent to the Feynman 
diagrammatic expansion for the logarithm of the determinant P, Q], although the zeta 
function approach provides a much simpler form of the final expression. Consider regulating 
the determinant with dimensional regularization. The perturbative expansion in powers of 
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the potential V is 

,21 \ (l\k+l 



( det [M + m 2 ] \ ^ (~l) fc+1 w 
Vdet [X free + m 2 ]/ ^ fc 

= 0-50+iO+- • (3-D 

where the dots denote insertions of the potential V. Alternatively, we can expand the 
dimensionally regulated determinant in partial waves as 

In ( det[A4 + m*] \ = f ( det [My + rrf] 

\<te[M*~ + m*\J £f \ det [.M^ + m 2 



i=0 

oo 



det 

^deg(Z;d) In /,(zm) . (3.2) 



In the Feynman diagrammatic approach |l8|, the first two order terms, and A^ 2 \ are 
separated out: 

ln ( det[i^ cc+ + n J 2 ] ) = £ de s( Z ;^{^/K im )-^/K^)]o(v)-[^/K^)]ocva) 



+A (D_Ia( 2 ) . (3.3) 

With these subtractions, the I sum is now finite, and the divergence lies in the dimensionally 
regulated Feynman diagrams and A^ . We now show how to relate the expression ()3.3|) 
to the zeta function approach. 

Using dimensional regularization, the first order Feynman diagram is 

A® = [ d d xV{x) \imG{x,y) . (3.4) 

J x^V 

The Helmholtz Green's function in d dimensions is 
We now use the Gegenbauer expansion |2 



K v {\x-y\) 



\x — yr 

1 y 1 ;=o 



2T(„) J2(l + u) Kl+ ;} r) CT (cos 9) , (3.6) 



where \x — y\ = ^/r 2 + (r') 2 — 2rr' cos6. As x — * y, noting that Cf^ 2 1 (1) 
we find 

(21 + d - 2) f l + d-3 

1=0 



l + d-3 
I 



drrV{r) J- 1 ( J_ 2) J f ^ J K l+d/2 ^(mr) I l+d/2 -i(mr) 
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= ^2deg(l;d) [hi fi(im)] 0(y) , (3.7) 
1=0 

in agreement with the 0(V) term in the iterative expansion for In fi(im) in (J2.25I) . 
Similarly, the second order Feynman diagram is 

A® = Jd d xJd d yV(x)G(x,y)V(y)G(y,x) . (3.8) 

Using the Gegenbauer expansion (J3.6|) . together with the identity 

r^(sin^) 2 ^r(cos^)^(cos^) = 5 IV , (3.9) 

Jo *! (« + ^) r M 

we find 

°° /•oo />r 

A™ = 2^deg(/;rf) / drr V{r)Kt +d/2 _ x {mr) / drVV(r')/f +(2/2 _ 1 (mr') 

oo 

= -2^deg(Z;d)[ln/f(im)] 0(va) . (3.10) 

1=0 

Thus, -|A( 2 ) agrees with the 0(V 2 ) term in the iterative expansion for \n fi{im) in (|2.25jl . 
when summed over I with the appropriate degeneracy factor. Therefore, the Feynman 
diagrammatic expression (|3.3|) is indeed equivalent to the zeta function expression (|2.13|L 
with dimensional regularization. 

Now _ aiso the hrate parts. In the Fe^ dia g _ ic approach Q, the 
finite renormalized logarithm of the determinant ratio is defined by the subtractions in 
()3.3|) . together with the finite renormalized form of the first two Feynman diagrams, 

A 0) 

and A^ 2 \ For definiteness we consider the d = 4 case, in order to compare with previous 
work Then the renormalized logarithm of the determinant ratio is 



111 ( det^e^] ) = £(' + {^W^) ~ MMW) - Miiim)] 0(ya) } 



+/1 fin 2 fin 



where we have used ()2.14|) to identify In fi(im) with ijj^(oo)/ipg^ e (oo). In [l8[, the subtrac- 
tion terms are defined as 

[\n fi{im)] 0{v) = /i[ 1} (oo) 
15 



[ln^(z'm)] 



o(y 2 ) 



(3.12) 



where h[ k \r) is the solution to the differential equation 



r ai 



dr' 



+ 



2ml' l+1 (mr) 
Ii +1 (mr) 



1\ d 
r J dr 



h\ k \r) = V(r)h\ K ^>(r) 



(fc-i). 



hf = 1 



(3.13) 



Comparing with (|2.25j) . we can alternatively express these subtracted terms as 



[ln/i(im)] 0(V Q = y drrV(r)K u (mr)I y (kr) 
o 



pn /i(im)] 



o(v 2 ) 



drrV(r)Kl(mr) / dr' r' V(r')I^(mr') 



(3.14) 



Finally, the finite contributions in (J3.ll)) from the first and second order Feynman diagrams 

|3] (note the small typo in equation (4.32) of 



in the MS scheme are 

,2 



m 



o 



A 



(2) 
fin 



1287T 4 



drr 3 V(r) 



dqq 3 V{q) 



^ a/ 4m 2 + q 2 ^ / \J 4m 2 + g 2 



y a/ 4m 2 + g 2 — q J 



(3.15) 



where V(q) is the four dimensional Fourier transform of the radial potential V(r). 

With the terms subtracted in ()3.11|) [evaluated using either (J3.13J) or (|3.14|) ]. the I sum is 
convergent. So this expression yields a finite answer for the logarithm of the determinant. 
On the other hand, the actual subtraction terms and counterterms in (|3.11j) are different 
from those in (jl.K)j) . even though the final net answer for the finite renormalized determi- 
nant is numerically the same. The difference between the two approaches is that in ()3.11|1 
one subtracts the full 0(V) and 0(V 2 ) dependence of In fi(im), given in (|2.25|) . and then 
compensates this subtraction with the Feynman diagram counter-terms whose finite part in 



the MS scheme are given in (|3.15|) . On the other hand, in the zeta function computation, 
we subtract just the asymptotic form of these first two Feynman diagrams, as in ()2.27j) . as 
is required to analytically continue the zeta function to s — 0. Subsequently, in the zeta 
function approach we notice that even this is an over-subtraction, as part of this asymptotic 
behavior cancels against C' as {0), leaving the final expression (|1.1U|) . The regularized form of 



1(3 



these zeta function subtractions is different from the regularized form of the Feynman dia- 
grammatic subtractions, but the associated counter-terms are also different, in such a way 
that the net result for the determinant is identical. This follows analytically from ([3.7)1 and 



so serves as an explanation of similar 

n n 

26l . |27| . However, we note that the 



([3.1U)1 , and can easily be confirmed numerically. This a 
effects noted in one-loop static energy computations 
zeta function expression ()1.10|) has a significantly simpler form, with the subtractions only 
requiring simple integrals involving V(r), while the subtractions in (|3.11j) require the more 
complicated integrals (|3.14|) [or, equivalently, solving the differential equations in ([3.13)1 ] . 
and also require the Fourier transform of the potential in ()3.15[1 . 

IV. CONCLUSIONS 

To conclude, we have derived simple new expressions, (ll.8j) - (jl.l()j) . for the determinant 
of a radially separable partial differential operator of the form — A + m 2 + V(r), generalizing 
the Gel'fand-Yaglom result (|1.3j) to higher dimensions. This greatly increases the class 
of differential operators for which the determinant can be computed simply and efficiently. 
Our derivation uses the zeta function definition of the determinant, but the same expressions 



can be found using the radial WKB approach of [16|,ll7|. Furthermore, we have shown how 
these expressions relate to the Feynman diagrammatic definition of the determinant based 
on dimensional regularization 0]. These superficially different expressions are in fact equal, 
although the zeta function expression is considerably simpler. 

A number of generalizations could be made. First, in certain quantum field theory appli- 
cations the determinant may have zero modes, and correspondingly one is actually interested 
in computing the determinant with these zero modes removed. Our method provides a sim- 
ple way to compute such determinants. For example, in the false vacuum decay problem 
arising in a self-interacting scalar field theory in <i-dimensional space-time, the 
pref actor for the semiclassical decay rate involves the functional determinant of the fluctua- 
tion operator for quantum fluctuations about a radial classical bounce solution $ c ;(r). This 
fluctuation operator is a radially separable operator of the form considered in this paper, 
but there is a cMold degenerate zero mode in the / = 1 sector, associated with translational 
invariance of the classical bounce solution. It is straightforward to generalize the analysis 
of 0] for the d = 4 case to other dimensions, to find that the net prefactor contribution 
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from these / = 1 zero modes (including the collective coordinate factor 3j, |32j) has a simple 
expression solely in terms of the asymptotic behavior of the bounce solution: 



S[^ci]\ d/2 Uet' [M {l=l) + m 2 ] 



-1/2 



2tt 



det 



JVl {l=\) 



(27T) 



d/2-1 



$ 1$", 

^oo \^cl 



d/2 



(4.1) 



Here the constant $oo is defined by the normalization of the asymptotic large r behavior of 



the bounce: <S>d(r) ~ $ no r 1-d / 2 K, 



d/2— 1 1 



Another important generalization is to include 



directly the matrix structure that arises from Dirac-like differential operators and from non- 
abelian gauge degrees of freedom. The Feynman diagrammatic approach is well developed 
for such separable problem s 13311 ; for example it has been applied to the fluctuations about the 
electroweak sphaleron ^is^. and to compute the metastability of the electroweak vacuum 
|3fil |. More recently, the angular momentum cut-off method has been used to compute the 
full mass dependence of the fermion determinant in a four dimensional Yang-Mills instanton 



background [I7l | . to compute the fermion determinant in a background instanton in the two 
dimensional chiral Hig gs m odel 3?j], and to address the fluctuation problem for false vacuum 
decay in curved space 38]. A unified zeta function analysis should be possible, as there is 
a straightforward generalization of the Gel'fand-Yaglom result (jl.3j) to systems of ordinary 
differential operators 
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